Collapse of Stiff Polyelectrolytes due to Counterion Fluctuations 
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The effective elasticity of highly charged stiff polyelectrolytes is studied in the presence of coun- 
terions, with and without added salt. The rigid polymer conformations may become unstable due 
to an effective attraction induced by counterion density fluctuations. Instabilities at the longest, 
or intermediate length scales may signal collapse to globule, or necklace states, respectively. In 
the presence of added-salt, a generalized electrostatic persistence length is obtained, which has a 
nontrivial dependence on the Debye screening length. 
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A polyelectrolyte (PE) is an ionic polymer which when 
dissolved in polar solvents dissociates into a long poly- 
mer chain (macroion), and small mobile counterions. Be- 
cause of the electrostatic repulsion of the uncompensated 
charges on the polymer, the chain is stretched out to rod- 
like conformations. Naively, one might expect that the 
more highly charged PEs are the more stiff. However, 
this tendency is opposed by the stronger attraction to 
the counterions, which may condense on a highly charged 
PE, giving it a much lower apparent charge |^,^|, and re- 
sulting in a lower stiffness. Nonetheless, in a mean- field 
(Poisson-Boltzmann) treatment of the counterions, the 
conformation of the PE is stretched. 

On the other hand, it is known experimentally that 
both highly charged flexible PEs (such as polystyrene 
sulphonate) |§g, and stiff PEs (such as DNA) §], can 
collapse in the presence of multivalent counterions to 
highly compact states. Rodlike PEs also form bundles 
under similar conditions To account for these effects, 
a relatively long-ranged attractive interaction capable of 
competing with the residual coulomb repulsion (of the PE 
backbones with compensated charge densities) is needed. 

To understand the origin of the attractive interaction, 
consider two like-charged substrates on which counteri- 
ons are condensed. As the two objects approach each 
other, the counterions may rearrange their positions. 
Any "correlated separation of charges" now leads to an 
attraction whose range is of the order of the "correla- 
tion hole" of the counterions on each substrate. Several 
different mechanisms may lead to correlated charge sep- 
aration: At low temperatures, the counterions form a 
Wigner crystal on the charged substrate, which leads to 
an attractive interaction with a range set by the lattice 
spacing ij^j. Another mechanism for charge separation 
is a specific binding of counterions to the substrate 
In this case the periodicity, and hence the range of the 
attraction, is dictated by the underlying structure of the 



substrate. Thermal fluctuations can also induce instanta- 
neous charge separations which inter-correlate on the two 
objects, leading to an attraction similar to the van der 
Waals interaction. Since dominant contributions come 
from fluctuations with wavelengths of the order of the 
separation of the objects, the range of the attraction is 
set by the distance between them; i.e. the interaction is 
long-ranged g|fi^. 

Recently, there have been a number of studies of the 
role of counterion condensation in the transition between 
extended and collapsed states of flexible PEs P, |l~3| . |i~I[ . 
The situation is more subtle for stiff PEs due to their 
intrinsic rigidity Jl5] [l7|] . Here we employ path integral 



methods [ [I2p8|j22( | to study the energy cost of deforming 
a stiff and highly charged PE in the presence of ther- 
mally fluctuating counterions. In particular, consider a 
chain of total charge N c and length L, with a micro- 
scopic persistence length £ p , and average separation a 
between charges on its backbone, in a neutralizing solu- 
tion of counterions of valence z. The PE is highly charged 
when a is less than the Bjerrum length £b = e 2 /eksT, 
where e is the dielectric constant of the solvent. (For 
water at room temperature, £b — 7.1A.) We calculate 
the effective free energy of a fluctuating PE as a pertur- 
bative expansion in the deformations around an average 
rodlike structure. The linear stability of this structure is 
controlled by a spectrum Z (k), as a function of the de- 
formation wavevector k. A negative value of £{k) signals 
an instability at the corresponding wavelength, leading 
to phase diagrams as in Figs. 2 and 3. In particular, in 
the absence of salt (Fig. 2), we find that counterion fluc- 
tuations cannot trigger collapse of a stiff PE, unless it 
has a microscopic persistence length less than a critical 
value of 



£p = A x a x z 



=r|x|l-4-]. 
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in which A is a numerical constant given later. For 
£ p < £p, there is a finite domain of intermediate lengths L 
of the PE for which a collapsed conformation is favored. 

For a low concentration of added salt, we find an effec- 
tive persistence length of 

j f , a(a/£Bf c 2 

P P+ 16zHl-a/zl B f(Kay^(±) 2*1* (i)' 

(2) 

where k~ 1 is the Debye screening length, related to the 
salt (number) density n via k 2 — Att£bti, and C2 is 
a numerical constant given below. The above expres- 
sion, which is a non-trivial generalization of the Odijk- 
Skolnick-Fixman (OSF) electrostatic persistence length 
jjjj, is a sum of repulsive (+) and attractive (-) electro- 
static contributions. It can therefore be negative under 
certain conditions, which we take as an indication of a 
conformational instability of the rodlike PE (tendency to 
collapse) . This occurs in the regions of the phase diagram 
indicated in Fig. 3, and only for persistence lengths less 
than a critical value given by an expression similar to 
Eq.(|l]). A similar 'softening' contribution to the rigidity 
of charged membranes, and associated conformational in- 



stability, is reported in Ref. |^0| . There are related insta- 
bilities caused by surface fluctuation-induced interactions 
between stiff polymers on membranes plf . 




FIG. 1. Transverse deformations of a rodlike PE. 

Without loss of generality, we orient the PE along the 
x-axis, and parametrize its configuration by a 3D vec- 
tor R(x). For small deformations (no overhangs) as in 
Fig. 1, we can set R(x) — (x,r(x)), where the 2D vec- 
tor r(x) measures transverse deviations from a rod for 
< x < L. For each configuration we would like to 
calculate a constrained partition function Z[R], by in- 
tegrating over all counterions positions in solution. As 
this is quite difficult, we instead generalize the approach 
of Ref. pH , and integrate over a fluctuating counterion 
charge density along the PE backbone, as |22| 



Z[R] s e - 0K "m = e - m f Vq{x) exp /_ I / dx (g(*) - go) 2 _ i_b f dx dx 

J | 2J o (8q) 2 2 J a a 



The first term, (3W[R\ = l p J Ax{d 2 x R) 2 /2, is the en- 
ergy cost of bending the PE, and j3 = l/fc^T. In the 
second term, the condensate charge is integrated over a 
gaussian distribution of mean q Q = a/z£s, and variance 
(5q) 2 — z(l — ajzl B) pd| , p2] |. A finite n anticipates the 
effect of added salt [[231. 

the effective Hamil- 



Using the methods of Refs. |12 



tonian can be calculated perturbatively in the deforma- 
tion field r(x). In terms of the Fourier modes r(fc), the 
effective free energy cost of deformations is 

0n e o ff [r(x)\ = \ J g £{k) |r(fc)| 2 + 0(r 4 ), (4) 



where in the absence of salt (k = 0) 



£{k) 



(qi/a 2 )£ B fc 2 ln(fc/fc ) 
[l + 2(8q) 2 (£ B /a) HL/a)} 2 



~B!k 3 



(5) 



(Despite its appearance, Eq.(^|) is not a Taylor expan- 
sion in k.) The rigidity term, £ p k 4: , ensures stability of 
the PE at large k (small length-scales). The first term 
describes the stiffening of the PE due to Coulomb inter- 
actions, which is most apparent at long wavelengths, and 
also involves a cutoff fco *C n/L pj| . Note that the fluc- 
tuations in charge reduce the strength of this term con- 
siderably, compared to the uniform case with (5q) 2 = 



p5| , p6| ). The new element is the term —Bik 3 which is 
due to the fluctuation-induced attractions [ p7t , and may 
lead to instabilities at intermediate k. The coefficient Bi 
can be obtained as an integral that depends weakly on 
{5q) 2 {£B / a), and that for typical values is well approxi- 
mated as B\ ~ ci/ln(L/et) 2 with c\ ~ 0.101. 

The rodlike configuration is unstable if any £(k) is neg- 
ative. The onset of instability is thus determined from 
£{k m ) — 0, at the minimum of the spectrum given by 
d£(k m )/dk = 0. In the limit of L 3> a, these con- 
ditions can be expressed as £ p = £ p , with the critical 
persistence length given by Eq.(]l|) above, with A s f ~ 
c 2 /ln 2 ( J L/a)ln[ci/2£ p /c ln 2 ( J L/a)]. For £ p < £ c pl there 
is domain of unstable modes for fc_ < k < fc + , where 

k± = [a/2£ p ln 2 (L/a)] (l±^Jl- £ p /£ c p ) . For finite 

chains, the values of k are restricted to k > ir/L, leading 
to the following three possibilities: (i) For k + < tt/L, the 
unstable modes cannot be accessed and the PE retains its 
extended form, (ii) For fc_ < tt/L < k + , the PE is unsta- 
ble at the longest scales. It will fold in a manner which 
we interpret as a precursor to collapse into a globular 
state, (iii) For tt/L < fc_, the PE is stable at the longest 
scales, but unstable at intermediate lengths. Although 
this analysis does not determine the fate of the unstable 
PE, we guess that this instability signals the formation of 
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a necklace structure p8| . These tentative identifications 
are depicted in the phase diagram of Fig. 2. The bound- 
ary of the collapsed phase is obtained from k± = tt/L, 



L/a 



a 7rln (L/a) 



ci 



(L/a)ln(n/Lk Q ) 
ATTz 4 (£ B /a) 3 (l-a/z£ B y 
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and has a nearly parabolic shape. 
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FIG. 2. Phase diagram for the salt-free case, for z = 4. 
The parameters a = 1.7 A, and fco = 10 _3 7r/L are used. 

For a low concentration of added salt (k ^ 0), the en- 
ergy cost of deformations is again given by Eq.(Q), but 
now the function £ (k) is analytic, with a well defined 
power expansion in k. The first term in the (properly 
regularized) expansion is Lpk 4 , from which we can define 
an effective persistence length given by 



L , 



4 [1 + 2(Sq) 2 (£ B /a) \n(l/ K a)} z (naf 



Eh 
Ik 



(7) 



Once again, L>i is obtained from an integral which has a 
very slow dependence on (5q) 2 (£ B / a), and which for typ- 
ical values, can be approximated by E>2 — C2/ln(l/«a) 
with C2 — 0.288. This leads to the expression for the per- 
sistence length in Eq.(|J). It is interesting to note that a 
similar result has been predicted for stiff polyampholytes 

The second term in Eq. (0) reproduces the OSF electro- 
static persistence length in the limit (Sq) 2 = |l9[] , with a 
reduced charge density go- Upon including the counterion 
fluctuations, (Sq) 2 ^ 0, there is a further reduction in 
this term. The final term is a negative contribution com- 
ing from the fluctuation-induced attractions. The latter 
reduces the 'rigidity' (effective persistence length) of the 
PE, which if negative results in conformational instability 
(collapse). Using this criterion, we obtain the phase dia- 



gram of Fig. 3. As in the salt-free case, the phase bound- 
ary has a nearly parabolic shape, with a maximum at the 
critical persistence length given in Eq. (Q) , with A as = c\ . 
Note that the maximal persistence length is much larger 
in the presence of salt, making it much easier to achieve 
conditions favorable to collapse. 
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FIG. 3. Phase diagram for the added-salt case, for 
different values of the counterion valence z. 

Rather than describing the formal steps of the method 
and approximations (which will be detailed elsewhere 
PH ) , we conclude with a qualitative summary of the na- 
ture of the results, and the range of their validity. From 
dimensional analysis, it is easy to show that unscreened 
Coulomb interactions make a contribution of £ B (k/a) 2 
to the rigidity spectrum £(k). In the PB solution, due 
to charge condensation the strength of this term is re- 
duced by a factor of q$ = (a/z£ B ) 2 . If the charge den- 
sity on the PE is allowed to fluctuate, it is further re- 
duced, and the Coulomb rigidity goes down by a factor of 
(Sq) 4 (£ B /a) 2 hi 2 (L/a), with (Sq) 2 = z(l - q ). However, 
these reductions do not change the overall sign which still 
prefers a rodlike structure. An attractive (destabilizing) 
contribution is generated by fluctuation-induced inter- 
actions, which are typically independent of microscopic 
parameters, and hence make a contribution of — fc 3 to 
£(k). Comparing to the leading Coulomb contribution, 
the latter corrections become important at short scales of 
order of a 2 /£ B . We thus may well question the applica- 
bility of continuum formulations to describe such a short- 
distance instability. In hindsight, the phase diagrams of 
Fig. 2 and Fig. 3 indicate that the prefactors involved in 
softening of the residual repulsion conspire to make the 
actual instability lengths quite large (~ z 4 £ B /a 2 ), and 
thus the continuum formulation should hold for a large 
portion of these phase diagrams |30|. Another poten- 
tial concern is the relevance of higher order loop correc- 
tions when the instability sets in. We are not able to 



3 



address this issue quantitatively, and it may be best to 
regard Eq.(||) as a variational approximation to the full 
problem. Finally, the instability analysis performed here 
only provides us with information concerning the onset 
of a conformational change. The final structure of the 
collapsed chain is naturally beyond this linear stability 
analysis pp3jj- 

It is known experimentally that different counterfoils 
with the same valence may behave differently as collaps- 
ing agents. The difference is usually attributed to modi- 
fications in the microscopic structure of the PE that take 
place upon binding of the counterions Q . It is thus plau- 
sible that the microscopic features that distinguish be- 
tween different counterions with the same electrostatic 
properties can be encoded in a single parameter, the 
microscopic persistence length, which normally depends 
only on the local microscopic structure of the PE back- 
bone. 
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